Desorption of n-alkanes from graphene: a van der Waals density functional study 
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A recent study of temperature programmed desorption (TPD) measurements of small n-alkanes 
(CjvH2jv+2) from C(0001) deposited on Pt(lll) shows a linear relationship of the desorption energy 
with increasing n-alkane chain length. We here present a van der Waals density functional study 
of the desorption barrier energy of the ten smallest n-alkanes (N — 1 to 10) from graphene. We 
find linear scaling with N, including a nonzero intercept with the energy axis, i.e., an offset at the 
extrapolation to N — 0. This calculated offset is quantitatively similar to the results of the TPD 
measurements. From further calculations of the polyethylene polymer we offer a suggestion for the 
origin of the offset. 

PACS numbers: 31.15.E-,71.15.Mb,71.15.Nc 



I. INTRODUCTION 

The increasing use of molecules on graphene and 
graphite surfaces for industrial applications calls for 
an improved atomic-scale understanding of the adsorp- 
tion/desorption structure and process. The n-alkanes 
are linear chains of hydrocarbons, short versions of 
the polyethylene (PE) polymer. Using temperature- 
programmed desorption (TPD) Tait et al. (Ref. [T]) mea- 
sured the desorption energy and desorption rate preex- 
ponential factor of n-alkanes on graphene deposited on a 
Pt(lll) substrate. The n-alkanes measured were short, 
with the number of C atoms N < 10. The desorption 
of n-alkanes from graphite surfaces was also measured 
by Paserba and Gellman [2H4] , and from various other 
surfaces by a number of other groups [5] , the surface ma- 
terials including metals (Ag, Au, Cu, Pt, Ru), oxides 
(AI2O3, MgO), and semiconductors (Si). 

In most of the alkane desorption measurements the des- 
orption energy was found to scale linearly with N for the 
short n-alkanes, but with a non-zero intercept with the 
axis of the desorption energy. The value found for this 
offset at N = was sometimes found to be unphysically 
large, several times larger than the scaling coefficient. 
In another study Tait et al. [6j [7] analyzed their own 
data for n-alkane on MgO(lOO) desorption. They allowed 
the desorption prefactor to vary with chain length and 
found the desorption energy offset to be non-vanishing 
but small, of the size of or smaller than the scaling 
coefficient. When the same group of authors analyzed 
their data of n-alkanes on Pt(lll) and on graphene (and 
also re-analyzed data from a number of studies by other 
groups for some of the above-mentioned surfaces) they 
found similar non-vanishing but small offsets for those 
desorption systems. 

In this paper we use the first-principles van der Waals 
(vdW) density-functional method 0(2], vdW-DF, to de- 
termine the n-alkane adsorption energy on graphene at 
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low coverages for short alkane chains (JV < 10). This 
adsorption energy can be compared with the experimen- 
tally determined desorption barrier energy values. As in 
the experimental studies in Refs. [TJ and [7] we find a 
close-to-linear growth in adsorption energy E a with chain 
length N, with a non- vanishing but small offset when ex- 
trapolated to N = 

E a = 7.23N + 6.44 [kJ/mol]. (1) 

Here and below we use the term adsorption energy (E a ) 
for the energy found in our theory calculations. This cor- 
responds to the desorption energy Eq of isolated alkane 
molecules on graphene. By Ed we denote the experimen- 
tal desorption energy, or desorption barrier, of an alkane 
molecule from partly covered graphene. In parts of the 
literature Ed is instead denoted Ai?j es . 

Contrary to analysis of the experiments, our calcula- 
tions of the adsorption/desorption energy do not involve 
an assessment of the desorption prefactor. Our values of 
E a are simply found from the differences in total energies 
of the system in the adsorbed and the desorbed states. 

The outline of the paper is as follows. In Section II we 
describe the method of computation, including conver- 
gence tests of sensitive computational parameters. Sec- 
tion III presents our results and discussions, including a 
discussion of the definition of a monolayer (ML) of cov- 
erage of n-alkanes on graphene. Section IV contains our 
summary. 

II. METHOD OF COMPUTATION 

The n-alkanes are linear, saturated hydrocarbon chains 
absent of branches, with the general formula CjvH 2 Ar+2, 
N > 0. Very long such chains (in principle infinitely 
long) are known as the PE polymer. In this paper we 
analyze the adsorption on graphene of the ten smallest 
n-alkanes {1 < N < 10), of H 2 , and of PE, all in the 
stretched form, which is the trans conformation. 

Our interest in the alkane desorption was sparked by 
the TPD experiments of Tait et al. p] and their analy- 
sis leading to the experimentally determined desorption 



energy. We determine the adsorption energy by use of 
first-principles density functional theory (DFT), employ- 
ing the method vdW-DF [5112] as detailed in several other 
publications [TUHT2"] but here with the vdW interaction 
treated fully selfconsistently [pj. We calculate the total 
energies of the adsorption system using the DFT pro- 
gram GPAW [H] with vdW-DF [SQ [§] in a Fast-Fourier- 
implementation |15j . 

Figure [T] illustrates the adsorbed n-pentane molecule 
on graphene, and the unit cell used in our calculations 
for N — 5. The lateral sizes of all the unit cells used are 
listed in Table |U 

For each of the adsorbed alkane molecules we deter- 
mine the optimal positions of the atoms by minimiz- 
ing the Hellmann-Feynman forces. These are derived 
from gradients in the self-consistently determined elec- 
tron density. This optimization also adjusts the in- 
tramolecular bond lengths to the most favorable value 
in the adsorption state. After this optimization we ob- 
tain the total energy of the adsorbate-graphene system, 

pvdW-DF 
near 

The GPAW code is an all-electron DFT code based on 
projector augmented waves [ST] (PAW) and using finite 
differences. In several of our previous vdW-DF applica- 
tions we used self-consistent calculations with the gener- 
alized gradient approximation (GGA) to determine the 
electron density and part of the total energy, followed 
by non-selfconsistent calculations to determine the to- 
tal energy within vdW-DF. This allowed us to focus on 
either the GGA need for accuracy in choice of computa- 
tional parameters and methods (for example the need of 
a fixed amount of vacuum in the total system [22H24"] ) 
or the vdW-DF need for accuracy in other parameters 
and methods (for example the need for fixating the lo- 
cal electron density grid [UJ [551 US])- When carrying 
out selfconsistent calculations for vdW-DF all of these 
requirements must be met in all calculations. Below the 
most important of these choices are described. 

We model the adsorption system by means of an or- 
thorhombic unit cell, periodically repeated in all direc- 
tions. The unit cell contains sufficient space in the 
plane of graphene for the repeated images of the n-alkane 
molecules to interact only little. As explained below, we 
make sure to explicitly subtract the small lateral interac- 
tion of the periodic images of the molecules from our re- 
sults [121 HZ] . In the direction perpendicular to graphene 
the introduction of ~ 19 A of vacuum above the adsorbed 
alkane ensures that no interaction across unit cell bound- 
aries takes place. 

For the smallest of the alkanes, methane (N = 1), we 
use a unit cell of size a g x 3a 9 x 23.0 A and for the 
largest alkane considered here, n-decane (N — 10), we 
use a fl 9 x4a 9 x 23.0 A unit cell. Here a g — V3a , 
with ao = 1.43 A, is the clean graphene lattice constant 
as found in our calculations by relaxing the lateral size 
of the unit cell. The unit cell sizes used for the other 
calculations (N = 2, ... ,9) are listed in Table |TJ 

We choose the real-space grid for representing the 
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FIG. 1: Illustration of n-pentane (N = 5) adsorbed on 
graphene. One unit cell is shown including a repetition of 
the graphene carbon atoms that are positioned on the unit 
cell boundary. C atoms are represented by large spheres and 
H atoms by small spheres. Figure created using XCrySDen 



wavefunctions in the PAW procedure to have a dis- 
tance less than 0.11 A between nearest-neighbor (nn) grid 
points. The (valence) electron density is represented on 
the same grid with additional grid points at half the nn 
distance, values obtained from interpolation of the elec- 
tron density and addition of compensation charges. The 
choice of these relatively dense (wavefunction and elec- 
tron density) grids is important for the quality of the 
evaluation of the nonlocal correlation contribution |25j . 

The Brillouin zone of the unit cell is sampled accord- 
ing to the Monkhorst-Pack scheme by means of a 2 x 2 x 1 
fc-point sampling. Increasing the fc-point sampling to 
4x4x1 changes E a less than 0.7 meV per molecule (0.07 
kJ/mol). We further make sure that the calculation is 
accurately converged with respect to total energies. We 
impose a convergence threshold such that the total en- 
ergy changes less than 0.1 meV per unit cell, or less than 
approximately 10 -6 eV per atom in the unit cell, in the 
last three iterations. This convergence threshold is sev- 
eral orders of magnitude smaller than the default settings 
of GPAW. 
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TABLE I: Adsorption (desorption) energies from theory (E a ) and experiment (Ed and Eq), center-of-mass distance from 
graphene d cm , area A of one alkane molecule in a full monolayer, unit cell used in calculations, and corresponding estimated 
coverages 8, for the small n-alkanes (CjvH2at+2, N = 1 — 10). In our calculations we use orthogonal unit cells and a graphite 
lattice vector a g = \/3 ao with ao = 1.43 A. The experimental values of Eq — corresponding to the limit of ML coverage and 
no defect sites — is found from Eq. |5| with the use of parameters given in Table IV of Ref. [T] The coverage for PE is found 
from an estimate of the PE-PE interaction distance in the PE crystal, as described in the text. 

This work Experiments 





N 


Unitcell 







E a 




A 


E d (0.5 ML) 


E 










[ML] 


[A] 


[kJ/mol] 


[eV] 


[A 2 ] 


[kJ/mol] 


[kJ/mol] 


H 2 





2^3 x 


3 




3.39 


6.8 


0.070 








methane 


1 


3^3 x 


3 


0.16 


3.64 


14.6 


0.152 


15 a 


14.1 


13.6 


ethane 


2 


3\/3 x 


3 


0.22 


3.80 


20.9 


0.216 


20.9 b 


24.6 


23.8 


propane 


3 


3\/3 x 


3 


0.28 


3.90 


27.7 


0.288 


27 a 


32.1 


30.6 


n— butane 


4 


3^3 x 


4 


0.26 


3.97 


34.6 


0.358 


32.7 C 


40.8 


38.9 


n— pentane 


5 


3\/3 x 


4 


0.31 


3.86 


42.8 


0.443 


39 a 






n— hexane 


6 


4x/3 x 


4 


0.26 


3.96 


49.6 


0.514 


45. 6 d , 44. 8 e 


63.0 


60.3 


n— heptane 


7 


4^3 x 


4 


0.30 


3.90 


57.7 


0.598 


51.6 / 






n— octane 


8 


5\/3 x 


4 


0.27 


3.89 


65.5 


0.679 


57.2 / , 57.7 d , 56.2 e 


72.6 


71.0 


n— nonane 


9 


5^3 x 


4 


0.30 


3.87 


73.0 


0.757 


63.5 / 






n— decane 


10 


5^3 x 


4 


0.32 


3.86 


80.3 


0.832 


69.0 / , 69.7 d , 68.9 e 


91.4 


84.5 


polyethylene 


(1) 5^3 x 


1 


0.21 


3.83 


7.2 s 


0.074 s 









"Linear interpolation of the experimental data available for other values of N, A(N) ~ 9 + 6iV A 2 . 

^Neutron diffraction data at submonolayer coverage, Ref. 1171 

c Neutron diffraction data at 11 K, Ref. fTHl 

d X-ray diffraction data at submonolayer coverage, Ref. 1191 

e Neutron diffraction data at submonolayer coverage, Ref. 1191 
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9 E a per C atom in PE. Each unit cell has two units of CH2. 



We determine the adsorption energy E a as the differ- 
ence in total energy between a system with an alkane 
adsorbed in the optimal geometry and a system with the 
alkane moved far away from graphene. It is well known 
that in the vacuum region of the system small spurious 
exchange energy contributions add to the total-energy 
[2"2Tf24] . To cancel these contributions we use the same 
unit cell for the calculation of the adsorbed state and 
for the desorbed state. We make sure that the height of 
the unit cell allows the fragments (graphene and alkane 
molecule) to be far apart within the unit cell. With unit 
cell height 23 A the maximum possible separation is ~ 11 
A, sufficient for the alkane to count as desorbed. 

The correlation energy E c in the total energy for the 
vdW-DF functional is split [28] into a nearly-local part 
E® and a part that includes the most nonlocal interac- 
tions Ef, 

E c = E° c + Ef . (2) 

In a homogeneous system the term E® is the correla- 
tion i^ DA obtained from the local density approximation 
(LDA), and in general [8] we approximate E® by E^ DA . 
The term E™ 1 vanishes for a homogeneous system. It 
describes the dispersion interaction. The form of E™ 1 is 



derived in Ref. [5] It is a truly nonlocal functional 

E?[n] = l[[ drdr'n(r)<Mr,r>(r') (3) 

given by a kernel </> which is explicitly stated in Ref. [51 

E^ 1 exhibits a small sensitivity to changes in the lo- 
cal real space grid, for example translations of the nu- 
clei positions by a non-integer number of real-space grid 
points [TTl [551 HU • The procedure to ensure high-quality 
results that include contributions from changes in the 
intra-molecular bonds is described in more detail in Refs. 
Upland |2T)1 and is reviewed here. 

We derive the adsorption energy in three steps. First 
we carry out a vdW-DF calculation of the adsorbed 
molecule, yielding the total energy E^^~ BF . Next we 
carry out a vdW-DF calculation of the molecule moved 
off graphene rigidly, which gives us the total energy 
-Ef£ " " . Finally, we let the molecule relax into the 
gas phase structure. From the last step we obtain val- 
ues of the total energy in the first (deformed for adsorp- 
tion) and in the last (gas phase) situation, E^^ lri and 
E^ E respectively. In this last relaxation study step we 
use the Perdew-Burke-Ernzerhof [3U] (PBE) variant of 
GGA. The adsorption energy E a , shown in Table [TJ is 
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then found from 

en _ pvdW-DF pvdW-DF 
— Ej a — £/ near — % ar 



E, 



PBE 
deform 



E, 



PBE 



(4) 



For the calculation of E^^~ DF we determine the 
optimal position of the alkane molecule adsorbed on 
graphene. The Hellmann-Feynman forces on the alkane 
atoms are minimized. All graphene atoms are fixed in 
space in this and all other calculations. 

For the calculation of £^ W ~ DF , we move the alkane 
molecule (rigidly) off of the graphene sheet by translating 
the molecule by 75 grid points, corresponding to adding 
~ 8.13 A to the graphene-alkane distance. This yields 
an alkane-graphene separation of approximately 11-12 A 
both between fragments within the unit cell and frag- 
ments in the vertically repeated images. By translat- 
ing the alkane an integer number of grid points the nu- 
clei locally maintain the same positions with respect to 
the grid, thus avoiding any possible effects of describing 
the electron density differently on a shifted grid [HO . 
This intermediate reference energy E^ W " DF is then sub- 
tracted from E^^~ BF . We use the same unit cell size, 
number of real space grid points, and number of /c-points 
as for the adsorption calculation. 

The use of an intermediate reference energy sub- 
tracts not only intra-molecular contributions (the alkane 
molecule is identical in the two calculations) but also any 
direct alkane-alkane interaction across unit cell bound- 
aries. This is usually a small contribution in our calcula- 
tions. 

The last set of calculations finds the energy gain 



pPBE 
^defori 



E P ® E of relaxing the isolated alkane molecule 



from the slightly deformed structure of the adsorbed 
state. This energy gain is not expected to have any long- 
range component. The E™ 1 term of vdW-DF depends 
slightly on the density grid points positions with respect 
to the molecule, and we therefore, as a matter of princi- 
ple, choose to perform this part of the calculation using 
instead the GGA variant PBE [3D]. 

It should be noted, though, that in this study it does 
not numerically matter to the energetics whether we use 
PBE or vdW-DF in this last relaxation-study step. 



III. RESULTS AND DISCUSSION 

The alkane molecules deform very little upon adsorp- 
tion, compared to their gas phase structure. In the gas 
phase, calculated with vdW-DF, we find for n-pentane 
the average C-C bond length 1.541 A. The C-C distance 
varies slightly along the carbon chain, with the smallest 
values towards the ends and the largest values around 
the center of the chain, but the difference only amounts 
to 0.001 A. For n-nonane the C-C bond lengths along 
the chain differ by 0.002 A, again with the largest val- 
ues around the center of the chain. We find similar bond 
lengths and bond length variations for the other alkanes. 

Values of the average C-C bond length found by ex- 
periment [3T] for n-alkanes with N = 2 to 7 are in the 



100 



90 



80 



70 



60 



50 



40 



30 



20 



10 



E d (0.5 ML), exper. □ 

E , exper. o 

E a , vdW-DF • 
PE & H 2 , 7.2 N + 6.8 




2 3 4 5 6 7 8 
Chain length N [number of C atoms] 



10 



FIG. 2: Desorption energy as a function of the length of the 
n-alkane chain. Solid points are our results, open points are 
from the TPD measurements by Tait et al. pQ. Linear regres- 
sion lines for the three sets of data points are also shown, in- 
cluding the extrapolations to N = 0. The dashed line has the 
slope from our calculations of the adsorption of PE, with the 
H2 adsorption energy added to represent the ends of the alka- 
nes. The shaded area encompasses the estimated (by Ref. H} 
errorbars of the 0.5 ML experimental data. 



range 1.526 A (n-propane) to 1.536 A (ethane). Our re- 
sults for the bond lengths thus deviate less than 1% from 
experiment. 1 

When an alkane molecule is adsorbed on graphene we 
find that the bond lengths arc only slightly affected. For 
pentane the adsorbed molecule (Fig. [I]) has an average 
C-C bond length 1.543 A, a change from the gas phase by 
0.002 A (or 0.1%), and the bond lengths are still larger 
towards the middle of the chain (1.544 A) compared to 
the bonds at the ends of the chain (1.542 A). For the 
other alkanes we find the average C-C bond length in the 
range 1.540 A (ethane) to 1.544 A (octane and nonane), 
and for all alkanes the bonds towards the middle of the 
chains are longer than those at the ends. 

Thus, structural changes caused by the adsorption are 
very small. Energetically, the changes are also very small: 
using PBE for the reasons stated in the previous section 
we find the difference in total energy E ( 



PBE _ 771PBE 
deform gas 



1 When comparing the values of the bond lengths with experi- 
ment it should be kept in mind that many exchange-correlation 
approximations, like the vdW-DF and also many of the GGA 
approximations, find covalent bond lengths that can deviate up 
to a few percent from the experimental values. 
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approximately 2 meV (0.2 kJ/mol) per C-C bond for all 
the alkanes studied here. 

All calculations of the adsorbed alkanes presented 
above are for an orientation with the alkane carbon skele- 
ton parallel to graphene. We started out the process of 
optimizing the atomic positions with the carbon skele- 
ton parallel to graphene, and all calculations reached 
an energetic minimum when parallel to graphene. To 
check if this is also the global minimum we further calcu- 
lated the total energy for alkane molecules that initially 
were oriented with their carbon skeleton perpendicular to 
graphene, i.e., rotated 90° around their axis. For these, 
we also find a local minimum, but with a total energy 
larger (less favorable) than for the configuration with the 
backbone parallel to graphene. For example for pentane 
we find that the loss of total energy going from the per- 
pendicular orientation to the parallel orientation is 48 
meV (4.6 kJ/mol). 

Table [I] lists the adsorption (desorption) energies ob- 
tained from theory by us and through TPD measure- 
ments by Tait et al. pQ . As shown in Figure [2] the cal- 
culated adsorption energy values grow linearly with the 
size of the alkane molecule, N, with an off-set comparable 
to that from experiments. Although the coverage of ad- 
sorbed alkanes in our calculation is 0.2-0.3 ML (with full 
coverage defined as described further below) the subtrac- 
tion procedure involving the two first terms in Q ensures 
that all direct alkane-alkane interactions are eliminated. 
Thus our results should be compared with experimental 
results for single alkane molecules desorbed from other- 
wise clean and defectless graphene (coverage ML), Eq. 

In Ref. [U the model used for describing the desorption 
energy Ed as a function of the coverage 9 and the number 
of alkane carbon atoms N is 



E d (6, N) = E (N) + y(N)6 + E dc{ (A)exp ( - 



7 dcf 



(AO 



(5) 

The Eq(N) is the contribution from a defectless surface 
(here: graphene) in the absence of adsorbate-adsorbate 
interactions. The term j(N)9 accounts for the increase 
in desorption energy due to the interaction with other ad- 
sorbates on the surface, and the third term describes the 
effect of defects in the surface. The model is introduced 
in Ref. 6]for n-butane on MgO(lOO). For general (small) 
?i-alkanes on graphite the parameters 7, -Edcf , & n d #def are 
given in Table IV of Ref. [TJ In our calculations graphene 
is defectless and there are no lateral interactions between 
molecules. Therefore our adsorption energies E a should 
be compared with the experimental quantity Eq, listed 
in Table U 

In Figure [2] are shown the experimental results Ed at 
6 = 0.5 ML, with error estimates within the shaded area 
(from Table III of Ref. [T]), and E at zero coverage, to- 
gether with our adsorption energies E a . Comparing E a 
with Eq we find that our theory adsorption energies devi- 
ate somewhat from the experimental results, with values 
from theory about 10% smaller than Eq. Hexane is an 
exception: our vdW-DF value E a is 18% smaller than 
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FIG. 3: Potential energy for n-pentane on graphene. The 
potential is calculated like — E a in Q but with the center of 
mass of pentane fixed at the distances d from graphene. The 
minimum of the curve corresponds to the adsorption energy 
E a and distance d cm . 



E . However, for hexane the experiment deviates from 
the linear growth with N whereas our theory result does 
not. 

The solid linear curves in Figure [2] are the linear re- 
gression curves for Ed, Eq, and E a . The experimental 
curves are described by Ed = 8.50A + 7.11 (Ref. [1]) and 
Eq = 7.96A + 7.46, and we find for the theory results 
the relationship E a = 7.23N + 6.44, all given in units of 
kJ/mol. 

The shaded area in Figure [2] shows the estimated errors 
in the values of Ed as provided by Ref. Q] it is reasonable 
to expect similar error estimates on the experimental Eq 
values but these are not available to us. 

In Figure [3] we show the potential energy curve for 
n-pentane as pentane is moved away from graphene, ob- 
tained with the vdW-DF functional. The points of the 
curve are obtained as described by Q, with three sets 
of calculations, however with the center of mass of the 
molecule kept fixed at the distance d above graphene. 
All internal atomic positions of pentane are allowed to 
relax. 

Figure [3] shows that the potential is shallow around the 
adsorption position. For example, moving the center of 
mass of pentane 0.1 A towards or away from graphene 
results in an energy increase of only about 1 kJ/mol, or 
10 meV. For our calculations we also report in Tablc|l]thc 
distance d cm between the center of mass of the molecule 
and the graphene sheet at the adsorption distance, which 
ranges from 3.64 A (methane) to 3.96 A (hexane). Be- 
cause of the shallow minimum there is some uncertaincy 
in determining d cm . 

Although we subtract the direct alkane-alkane interac- 
tions from the adsorption energy it is in principle possi- 
ble that small indirect alkane-alkane interactions remain 
if the molecules are placed too close, e.g., closer than in 
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a monolayer coverage. Such indirect interactions would 
affect the atomic positions of alkane from nearby alka- 
nes. We must therefore make sure that the alkanes are 
sufficiently far apart that such indirect interactions are 
negligible. To quantify this, we determined the coverage 
of alkanes in our calculations in units of ML, for each of 
the alkanes as well as for H 2 and PE. 

We describe the coverage of molecules on graphene by 
fractions 9 of a molecular ML. One ML is a one-molecule 
thick coating of a surface, as found by experiments. In 
order to find 9 for our calculations we need to know for a 
full ML how many molecules cover a specified area of 
graphene, or how large an area A does one molecule 
cover, in average. Couto et al. |32j found by means of 
STM that for various n-alkanes adsorbed on graphite the 
coating layer is highly ordered. The ordering at such 
high coverage is affected both by the adsorbate-adsorbate 
and the adsorbate-substrate interactions. Disordered ar- 
rangement is only activated above a critical temperature. 
For most of the molecules we find the definition of 1 ML 
from experiments reported in the literature |17H20j . but 
for methane, propane, and pentane we use an estimate 
based on a linear interpolation of the experimental data 
available for other values of N, A(N) 9 + 67V A 2 . From 
the values of A(N) and the sizes of our unit cells we cal- 
culate the coverages 9(N) used in our calculations. 

For PE we estimate 9 from the PE-PE interaction dis- 
tance in the PE crystal. In an earlier study [53] one of 
us identified the optimal PE crystal structure within a 
vdW-DF characterization. From that study we estimate 
the optimal centerline-to-centerline distance between the 
polymers to be 4.5 A. The separation of the PE polymers 
used in the present study is 5\/3a g 21 A. This gives a 
coverage for PE 6> PE « 4.5 A/21 A « 0.21. 

As summarized in Table ||J the coverages used in our 
study are always less than a third of a ML, with values 
9 = 0.16 — 0.32, sufficiently sparse for indirect interaction 
effects to be neglected. 

Almost all the available experimental results on alkane 
desorption from various surfaces derive from TPD mea- 
surements measuring the desorption rate r. In order to 
extract the desorption energy from r the preexponen- 
tial desorption rate v was earlier often assumed to have 
the value 10 13 s _1 . This value is accepted as a reason- 
able value for first-order processes in surface physics of 
atoms and is derived from traditional transition state the- 
ory. However, the more complex processes of molecular 
desorption are not necessarily as well described by that 
particular value, nor more generally by a value that is 
constant for all n-alkanes. 

The desorption rate r may be described by the Polanyi- 
Wigner equation 

r(0,T) = -j-(6,T) = i/(0, T)9 n e~ Ed{e)/kBT (6) 

for 7ith order desorption, here n = 1. Assuming a con- 
stant value of v for the small (N < 12) n-alkanes the 
TPD desorption rates give linear growth in Ed with num- 



ber of alkane segments N but with a very large offset [5J 
at N — 0. The offset is much larger than the segmen- 
tal increment in Ed- Lei et al. speculated |34| . and Tait 
et al. showed from analysis of TPD experiments [TJ [7] , 
that v takes other and varying values in alkane desorp- 
tion. This was shown for various surfaces like graphite, 
Pt(lll) and MgO(lOO). By treating v as a fitting pa- 
rameter along with Ed, modified and varying values of v 
are found. Such analysis leads to a more modest value of 
the offset of Ed at N = 0, at the size of or smaller than 
the segmental increment in Ed [JJ [7] . 

In particular, it was found [1 that on graphite v varies 
from 10 130 s" 1 for methane to 10 17 8 s _1 for n-decane. 
Thus the small molecules have a desorption prefactor 
similar to that from theory for atoms, whereas the pref- 
actors for the larger molecules deviate strongly from this. 
Taking these variations into account the desorption en- 
ergy offset at N = is reduced to 7.11 kJ/mol, with 
a segmental increment in Ed of 8.50 kJ/mol. Similar re- 
sults were obtained for small- N n-alkane desorption from 
Pt(lll) and MgO(100). 

It thus seems that the previously published large values 
of the N = offsets can mostly be explained [H [51 E] as 
an effect of not allowing v to vary for the small n-alkanes. 
Nevertheless, an offset of a smaller size does remain even 
in the re-analyzed data. 

In the literature the origin of the offset has been de- 
bated [TJ [5J 153] . Even though the values of the offsets 
may be reduced as discussed above, also the remaining 
offset begs an explanation. Lei et al. summarize the dis- 
cussions by listing a number of suggested reasons: (i) the 
different binding to the surface of the methyl end groups 
(-CH3) compared to the methylene segments (-CH2-); (ii) 
the effect from the chain length dependence on the polar- 
izability of the alkanes; (iii) the effect of needing different 
temperatures for the various alkanes for measuring r; (iv) 
possibly the desorption process cannot be described as a 
first-order process, e.g., if the alkanes adsorb in islands 
or other structures; (v) possible lattice mismatch of the 
alkanes with the surface; and finally (vi) chain length de- 
pendence of v. The latter suggestion reduces the offset 
to a more modest value, as discussed above. 

Without going into details of all of the above- 
mentioned suggestions we note that our calculations are 
in a sense more direct than the desorption energies de- 
rived from the TPD measurements. In our calculations 
the preexponential factor v is not involved, temperature 
variation is not an issue, and we do not let the alkanes 
adsorb in islands. Our results are in agreement with the 
results presented in Ref. UJ where the approach of a vari- 
able desorption prefactor was used. In particular, our 
theoretically calculated value of the offset agrees very 
well with that obtained by Tait et al. Here we present a 
simple model study to discuss the suggestion (i) of end- 
group effects. 

Our calculated adsorption energy for PE corresponds 
very well with our similarly calculated adsorption energy 
per segment of the alkanes (when neglecting the offset). 
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PE is similar to the alkanes, but it does not (at least 
not ideally) include methyl end groups. In our calcula- 
tions we describe PE adsorbed on graphene by period- 
ically repeating two CH2 units, thus explicitly avoiding 
end groups. We find (Table |l| the adsorption energy per 
methylene unit in PE, 7.2 kJ/mol, which corresponds 
very well to the energy 7.4 kJ/mol we find per (methy- 
lene or methyl) unit for small n-alkanes. 

It is natural to expect that the two extra H atoms at- 
tached to the ends of the alkane molecules (in the methyl 
groups) also contribute to the adhesion, thus affecting 
the offset in the adsorption energy. We present a calcu- 
lation of a H2 molecule adsorbed on graphene to test a 
hypothesis of simple additivity of N methylene segments 
(-CH 2 -) and two additional H atoms. This is thus an 
even simpler model for n-alkane than adding methyl to 
the ends of a string of methylene segments. 

Our calculation of H2 on graphene yields the adsorp- 
tion energy 6.8 kJ/mol. In the calculated curve for E a 
the offset is 6.44 kJ/mol. Our results for PE and H 2 fit 
nicely to this simple additivity model. The curve with 
slope derived from PE adsorption and offset derived from 
H2 adsorption is plotted in Figure [2] (dashed line) . 

Arguments raised against the end-group explanation 
have been that experiments [34] for cyclic alkanes on 
Cu(III) and Pt(Ill) also show an offset for extrapola- 
tion to N even though the cyclic alkanes do not have any 
end groups. However, those results were extracted using 
fixed values of v and yield large offsets (36 kJ/mol for Pt, 
19 kJ/mol for Cu) both for the cyclic alkanes and their 
linear equivalents. In that analysis the effect directly on 
the desorption barrier from the end groups is estimated 
to 2 kJ/mol per linear alkane for both Pt and Cu surfaces. 
We cannot judge whether the value of the cyclic-to- linear 
difference in desorption barriers [34j , 2 kJ/ mol, would re- 
main after a re-analysis of the desorption energies with 
more variation of v, along the lines of those of Tait et al. 



IV. SUMMARY 

We present a computational study of the adsorption 
of small n-alkanes on graphene using the van der Waals 
density functional method vdW-DF. Recent desorption 
experiments pQ have shown desorption barriers growing 
linearly with the size of the alkane molecule, but with 
an offset in the limit of zero length. Here we reproduce 
in our calculations the linear dependence on the alkane 
length including an offset the same size as obtained by 
the experiments. With the help of our calculated adsorp- 
tion energy of polyethylene and H2 molecules we argue 
that a simple additivity assumption of alkane methylene 
(-CH2-) units plus two extra H atoms for the alkane ends 
explains the size and origin of the energy offset very well. 
Summing up, our calculations thus give support to the 
suggestion that the offset measured in n-alkane desorp- 
tion experiments (after correction for effects of varying 
v) can be explained by the n-alkane end groups being 
different from the methylene segments of the n-alkanes. 
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